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1. Introduction 

The idea of achieving enantiomer separation on the basis of specific transport properties 
of chiral objects in flowing fluids [f, 2, 3] is very attractive because it opens the 
perspective to tackle this also industrially relevant task by purely physical means without 
the need of specific chemical agents. 

Howard, Lightfoot and Hirschfelder [4] were the first to suggest the use of 
transport properties of asymmetric objects in order to separate differently crystallized 
enantiomers. Later, de Gennes discussed a mechanical way of separating chiral crystals 
[5]. Only recently, a lateral drift of screw-like objects in shear flows predicted by Brenner 
[6] was experimentally demonstrated for bodies with a linear extension of the order of 
centimeters [7], millimeters [8] and micrometers [9]. The direction of this drift, resulting 
from the coupling of the rotational and translational motion of a chiral object in a fluid, 
differs for objects with opposite chirality. The mutual interactions between flowing fluids 
and immersed chiral objects has been the subject of various recent experimental studies. 
The generation of flow patterns by actively rotating chiral objects was investigated in 
Ref. [10]. The control of actively moving nano-structured propellers was studied in 
Ref. [11]. In Ref. [12] the formation of chiral assemblies made of achiral components in 
a vortex flow is described. 

Theoretical studies of the influence of the chirality of a macroscopic object on 
its transport properties were performed for sedimentation [13] and for shear flows 
[14]. Chiral planar three-atomic molecules were studied in clover leave vortex flows 
[1] and in flows in channels of various geometries [15, 16]. In these two-dimensional 
studies fluctuating forces were taken into account in order to model the action of 
thermal fluctuations exerted on the molecules by the surrounding fluid. All studies 
[1, 13, 14, 15, 16] though are restricted to rigid objects. 

Watari and Larson considered three-dimensional, four-atomic molecules which, in 
equilibrium, assume the form of a regular tetrahedron, which as such is achiral with 
respect to geometry [17]. When put into a shear flow these molecules are distorted 
in an asymmetric way due to differently deformable bonds between the atoms. As a 
consequence these molecules experience a lateral drift. 

In the present work we study the motion of four-atomic molecules with different 
equilibrium bond-lengths but identical bond-strengths, see in eq. (2) below, being 
suspended in a uniform shear flow. The main mechanisms leading to a lateral drift of 
the molecules, i.e. to a motion perpendicular to the flow direction and to the direction 
in which the flow speed changes, can be summarized as follows. The shear field of the 
flow leads to a rotation of the molecule which then, like a small propeller, moves in, or 
opposite to the direction of the vorticity of the flow. We demonstrate that the rotation 
is essentially determined by the advection of the atoms with the fluid motion and only 
little influenced by hydrodynamic interactions between the atoms. In a uniform shear 
flow a coupling between rotational and translational motion, however, only is possible 
due to hydrodynamic interactions between the atoms. 
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As for the motion of molecular motors [2], the direction of the lateral drift is 
difficult to predict. In general it depends on all parameters of the system and also 
on the initial conditions with which the molecule is started. These initial conditions 
may belong to different domains of attraction, eventually leading to different attractors, 
which may display different forms of rotation of the molecule, and consequently different 
drift behaviors, even with opposite signs. 

Due to the deformability of the considered molecules, chiral attractors may exist 
for molecules with achiral equilibrium configurations and consequently also achiral 
molecules may experience a lateral drift displaying a dynamical symmetry breaking. 
The lateral drift, however, vanishes on average because for a molecule with achiral 
equilibrium configuration chiral attractors can only come in mutually mirror symmetric 
pairs. 

The paper is organized as follows. In Section 2 we introduce a class of "spiral" 
molecules, the geometry of which is described by only three parameters, and specify 
their equations of motion in terms of bond strengths and hydrodynamic interactions. 
Further, we introduce the most important quantities needed for the analysis of the 
numerical solutions of the equations of motion (Section 2.2). In Section 3 we present 
the results concerning the center of mass motion (Section 3.1), the rotational motion 
(Section 3.2) and the internal motion (Section 3.3) with particular emphasis on the 
dynamical symmetry breaking (Section 3.3.1). The role of hydrodynamic interactions is 
discussed in Section 4. The paper closes with a summary and an outlook in Section 5. 

2. The model 

2.1. Spiral molecules and equation of motion 

In the present case study we consider "molecules" consisting of four identical "atoms" 
modeled by spheres of radius a. The equilibrium configuration of a free molecule, i.e. 
a stable configuration with balanced bonding forces between the atoms, is assumed to 
be of helical form. This form is constructed by positioning the centers of the spheres 
on a cylinder mantle of radius p. Given the location of the first atom, R\, the next 
one, i?2, is found upon rotation by an angle tp about the cylinder axis combined with 
a translation by h parallel to the cylinder axis, and so on with the third and forth 
atom, yielding the respective positions -R3 and -R4. Hence, if the cylinder axis coincides 
with the z-direction of a coordinate system and the first atom lies on the x-axis of the 
coordinate system, then the following vectors point to the centers of the atoms sitting 
on a spiral: 

R n = pe x cos[(n — l)ip\ + pe y sin[(n — l)ip\ + h(n — l)e z , n — 1, 2, 3, 4.(1) 

We model the bonding forces between the atoms of such a "spiral" molecule by 
nonlinear, so-called FENE-Fraenkel (Finitely Extensible Nonlinear Elastic) springs [18]. 
If r is the vector connecting two atoms, then the force acting between these atoms is 
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given by 



f = k(r — I) 



(r + -l)(l 



r_) r 



for r_ < r < r 



+ > 



(2) 



(r+ — r)(r — r_ / 

where r = \r\ is the distance between the atoms, / the equilibrium length of the spring, 
k the spring-constant describing the Hookean behavior of the spring for small extensions 
about the equilibrium length, and s < 1 is a nonlinearity parameter of the spring. At 
the extensions 

r± = (1 ± s)l (3) 

the force diverges so that the bond-length between two atoms is restricted to the 
indicated, finite range. 

In contrast to Ref. [17], we choose identical force parameters k and s for the bonds, 
only the equilibrium bond lengths I may differ. Hence, a possible chirality of the 
equilibrium configuration of a molecule is solely determined by its geometry: According 
to the standard definition, a molecular configuration is achiral if a movement of the 
molecule, i.e. a combination of translations and rotations in three-dimensional space, 
exists such that it matches its mirror image; in any other case the configuration is chiral. 

When the molecule is suspended in a flowing fluid, advective forces act on the 
individual atoms, which in addition mutually influence their motions by hydrodynamic 
interactions. The total action of the fluid on the atoms will be described by a mobility 
tensor H n;rn relating the forces acting on the mth particle to the velocity of the nth 
particle. The deterministic, overdamped motion of the molecule in a fluid flow with the 
velocity field v(r) is then given by 



r n = v(r n ) + ^2 Hn >™ ' f" 



(4) 



m=l 



where f m is the total force acting on the mth atom, i.e. the sum of the FENE-Fraenkel 
forces (2) exerted by the other three atoms, and the dot denotes differentiation with 
respect to time. For sufficiently small spheres representing the atoms, inertial terms 
can be neglected [19]. We shall describe the mobility by the Rotne-Prager tensor which 
includes the hydrodynamical interactions up to second order in the ratio a/r n m of the 
atomic radius a and the distance r n>m between the nth and the mth atom. It is given 
by [20] 
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where 1 denotes the unit tensor in the three-dimensional space, 7] the viscosity of the 
fluid, and r n m = r m — r n the vector pointing from the position of the nth to that of 
the mth atom. 
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Figure 1. Uniform shear flow: The velocity is independent of x, points in z direction, 
and varies its magnitude proportionally to y. 



In passing we mention that thermal noise leads to additional random contributions 
in the equations of motion, see Ref. [21], yielding 

4 

V n = V(r n ) + H n,mfm + 2k B T <Tn,m " £ m (0 , (7) 
m=l 

where £ n (t) are mutually independent three-dimensional Gaussian, white noise vectors, 
satisfying 

(Ut)Zm(s)) = 5n, m 15(t - S) . (8) 

The fluctuation-dissipation theorem requires that the coupling to the noise, <7jj, is 
related to the Rotne-Prager tensor according to 

4 

Hn,m ^ ^ ""n,fc ' j,m ■ (9) 

k=l 

In the present work though we neglect all influences stemming from thermal noise, 
and, hence restrict ourselves to the discussion of the deterministic equations of motion 
in (4). We shall consider a uniform shear flow given by 

v x = v y = 0, v z = ny . (10) 

where k denotes the shear rate, see Figure 1. Even in this seemingly simple case 
the deterministic equations of motion can only be solved by numerical means. For 
that purpose we employed a LSODA code of variable order and step-length [22, 23] as 
available from the open source SciPy package [24]. 

In the numerical treatment we used dimensionless variables where lengths are 
measured in units of the atomic radius a and times in units of the inverse shear rate k^ 1 . 
Apart from six parameters which fix the equilibrium geometry of a molecule consisting 
of four atoms in general, or three in the present case of a spiral molecule as defined by 
eq. (1), two further parameters are needed for the characterization of the forces. As 
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such we will choose the nonlinearity parameter s and the dimensionless spring constant 
kf = k I '(6irr)a 2 k) . With this choice the deterministic equations of motion read 

4 

•En Sx n -\- kf ^ ^ Hn,m ' F m , (H) 
m=l 

where S is the shear tensor, which has vanishing matrix elements apart from the z-y 
element, being unity, and F m denotes the dimensionless force exerted on the j-th particle 
by the adjacent particles, reading 

, J- V J/ n,m/ ^n,m) Jj n > m 

for (1 s)^n,in *^ %n;rn *^ (1 "I - S)(n,m ■ 

Here T* n /o, x n ^ m 33 n 3! m) %n,m |^?i,m| and ( n ;m |-^n ■^m|/^ denote 

the dimensionless positions, distance vectors, the actual distances and the equilibrium 
distances between atoms n and m, respectively. 

Possible initial conditions for the integration of the equation of motion (11) are 
given by any set of atomic positions that obey the restrictions imposed by the FENE- 
Fraenkel forces (13). Here we confine ourselves to equilibrium configurations of the free 
molecules. For further details of how the orientations are chosen, see Sect. 3. 

2.2. Characteristics of the mechanical molecular states 

In the framework of the overdamped dynamics of eq. (11) the mechanical state of the 
molecule is uniquely described by the positions x n (t) of the four atoms. A reduced 
description is provided by the center of mass of the molecule consisting of four atoms 
of equal mass, consequently reading 



X 



z5>- ( 13 ) 



The spatial orientation of the molecule can conveniently be characterized by the three 
principal axes of the symmetric tensor of second moments of the atomic positions, M, 
defined by its elements 



'',.7 



ij>? (14) 



where denotes the i-th Cartesian position component of the n-th atom. The 
orientations of the principal axes are determined by the normalized eigenvectors Ui and 
their magnitudes by the corresponding eigenvalues \ of the tensor of second moments. 
Hence, they result from the eigenvalue problem 

Mui = XiUi , (15) 

where M = (M id ). 

The velocity vector u> of the instantaneous rotation of the molecule in the frame co- 
moving with the molecule's center of mass is determined by the orientation and velocity 
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of any body-fixed tripod. Here we use the principal axes of the second moments to 
obtain 



= -^-^2u i xu i , (16) 



where the dot denotes a derivative with respect to time and u x v the vector product 
of the vectors u and v. By means of the eigenvalue equation (15) one may express ttj 
in terms of the time derivative of M, to yield 

- 2 E^ x ^ ( 17 ) 

where u ■ v denotes the scalar product of the vectors u and v. 

As a measure of the deformation of a molecule we consider the ratio of volumes 
taken in the distorted state and the equilibrium shape. As molecular volume one may 
assign the square root of the product of the magnitudes of the principal axes, i.e., 

V = (AxA 2 A 3 ) 1/2 . (18) 

Whether a molecule is chiral cannot be decided on the basis of its second moments 
of positions. For this purpose, we use the isotropic chirality index Go, which is invariant 
under translations and rotations of the molecule [25], in its dilation invariant form. For 
a molecule consisting of four atoms of equal mass and size it becomes 

Go = 1 ^ XiJ ' ( Xk > 1 X ( Xi ' j ' Xj ' k ^ ( Xj ' k ' Xk ^ (19) 

The sum is extended over all sets i,j,k,l, taken from the permutations of {1,2,3,4}. 
This index has the general property that it assumes the same absolute value but opposite 
signs for enantiomers, i.e. pairs of configurations that are mirror images of each other. 
Hence, achiral configurations have a vanishing chirality index. 

For a spiral molecule, the chirality index is a function of the relative rise p/h and the 
screw angle ip. It is an antisymmetric function about the line p = ir, see Fig. 2(a), and 
possesses a singularity at <p — 2n/3 and h/p = 0. Along the lines h/p pa — (p — 2ii/3) for 
ip < 27r/3, and h/p pa 0.95(</? — 27r/3) for <p > 27r/3, the chirality index forms a ridge and 
a valley, respectively. Approaching the singularity along the ridge, the chirality index 
converges to the approximate value 0.289 while it approaches the opposite value —0.289 
along the bottom of the valley. 

Since all configurations with h/p = are planar, and, hence achiral, the chirality 
index exactly vanishes on the line h/p = 0. In the remaining part of the y?-/i//)-plane 
the chirality index has three extrema on each side of the symmetry line ip = tt. Local 
extrema are located at ip pa 0.8908tt, h/p pa 0.2073 and ip pa 0.9291vr, h/p pa 0.584. 
The absolute maximum is at ip ~ 0.6567T, h/p pa 2.516, see Fig. 2(a). The chirality 
index vanishes on the solid and dashed curves in Fig. 2(b). These curves result from 
constraints on the bond lengths (see Appendix). Along the dashed curve emerging from 
(p = 7r/2, h/p = 0, and ending in p = n, h/p = 0, the configurations possess a group of 
four and another one of two equally long edges, see Fig. 3(b). Since the corresponding 
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configurations possess a symmetry plane, they are achiral. The faces of these tetrahedra 
are made of identical isosceles triangles. Also the configurations on the curve starting 
at tp — 27r/3, h/p = and reaching tp — tt, h/p = 2/V3 are achiral; the faces of these 
configurations consist in two identical equilateral and two identical isoscele triangles; 
consequently five of the six edges have the same length, see Fig. 3(c). The analytic 
forms of the achiral curves are derived in the Appendix. Examples of different types 
of configurations are displayed in Fig. 3. There exists a third node line of G which, 
however, does not correspond to achiral configuations. We could not find an analytic 
expression for this line which emerges from tp « 0.467T, h/p — and ends in p 0.88-7T, 
h/p = 0. The three node lines of the chirality index meet in a single point which 
corresponds to a regular tetrahedron. For the fact that a vanishing value of the chirality 
index does not necessarily correspond to an achiral configuration, see in Ref. [26]. 

3. Results 

For our investigations, we selected spiral molecule configurations on four cuts in the 
parameter plane displayed in Fig. 2(a). Three of these cuts (A, B, C) are taken 
at constant relative rise h/p = 0.2073,0.584, and 2.5162, respectively, each one for 
p G (0, 7r). For the two smaller values of h/p we chose a relatively large radius p = 7 
in order to avoid any collision of the atoms within a molecule, and the smaller radius 
p = 2 for the large relative rise. The vertical cut D is determined by the screw angle 
ip = 27r/3, radius p = 3 and relative rise h/p G (0, 5). 

As initial configurations we took the equilibrium configurations according to eq. (1), 
rigidly translated such that the center of mass falls into the origin. Finally, rigid 
rotations of the molecule about its center of mass were performed. Between 500 and 
5000 rotations were sampled from the uniform distribution on the three-dimensional 
rotation group [27], i.e. from the Haar measure of SO(3). We typically let run each 
trajectory for 4000 time units [/c -1 ]. In the majority of the cases the trajectories have 
relaxed toward a stationary regime long before the time span of 4000. In this state 
the center of mass acquires a finite velocity in ^-direction, depending on the precise 
initial position of the molecule, and the y-component of the velocity may perform small 
oscillations but vanishes on average. The lateral motion in the direction of the vorticity 
coinciding with the x-direction in the present case is characterized by an average velocity 
superimposed by bounded modulations which were always periodic. Fig. 4 displays the 
components of the center of mass for the same molecule and its mirror image in four 
different asymptotic motional states. 

3.1. Center of mass motion 

The fact that the rc-component of the center of mass moves at a finite average velocity is 
apparently related to the chirality of the considered molecule: For any achiral molecule 
in the considered shear flow the opposite x-directions are equivalent. This symmetry is 
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Figure 2. The chirality index Go defined in eq. (19) is displayed for spiral molecules 
as a function of the screw angle ip and the relative rise h/ ' p, see eq. (1). On the panel 
(a) the symmetry Go — > —Go upon a reflection on the line tp = n is clearly visible. The 
four black lines labeled by A - D indicate four cuts through the parameter plane, see 
the text in Section 3 as well as Fig. 5 below. The lines A and B are defined by local 
maxima of the chirality index. On C and D the absolute maximum of Go is located. 
The panel (b) represents a blow-up of the region marked in the panel (a) by the white 
lines. Note that for the purpose of better visibility different color scales are used in 
the two panels. The black lines in (b) indicate those curves on which two of the three 
edge lengths l(k), k = 1,2,3 between atoms with positions R n and R n +k are equal, 
see the Appendix. The solid line represents achiral configurations with 1(1) = 1(2) 
which have five equally long edges. It starts at the singularity of Go at tp — 2w/3 
and h/p = and ends at <p = An/ 3 and h/p — 0. Along the dashed line 2(1) = 1(3) 
holds. These configurations have two and four equally long edges and are also achiral. 
On the dashed-dotted line, where 1(2) = 1(3), the configurations possess two triples 
of equally long edges, and the chirality index is non-zero everywhere except at the 
points with h/p = and at the crossing point of the three lines which corresponds to 
a regular tetrahedron. Note that there is still a third curve recognizable by the green 
color on which the chirality index vanishes in spite of the chirality of the corresponding 
configurations. 

broken by the presence of a chiral molecule and, according to the Curie principle [29], the 
molecule will perform a directed motion in either of the two directions. This argument is 
further corroborated if one launches the chiral partner of a considered molecule into the 
shear flow with initial atomic positions obtained from the original ones by mirroring on 
the y-z plane. Then one indeed finds the molecule moving in the opposite x-direction 
as illustrated by the panel pairs ((a), (ma)), ((b), (mb)), etc. in Fig. 4 displaying the 
corresponding pairs of mirror images. 

Inserting a molecule into the shear flow with randomly chosen orientations as 
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Figure 3. Panel (a) illustrates a generic spiral molecule with three and two edges of 
equal lengths and a sixth edge of different length. The edges of equal lengths are drawn 
correspondingly. Panel (b) presents a configuration on the dashed curve displayed in 
Fig. 2(b) with two and four equally long edges while the configuration in panel (c) 
has five equally long edges representing a configuration on the solid curve depicted in 
Fig. 2(b). The right part of the figure represents the corresponding mirror images. 
The cases (b) and (c) are achiral. 

described above, one typically finds between one and four different motional states 
of the center of mass, which are asymptotically reached at large times. We determined 
the corresponding average velocities (v x ) by a least-square fit of the x-components of 
the center of mass to a uniform motion, yielding 

N 

12 N 

(fx) = T77T7 ^ttt )X x (t + nr). (20) 

N 1 N(N + l)(JV + 2)r ^ 2 ' y ' y ' 

In most cases we used t = 3000, iV = 1000 and r = 1. Since the different motional 
states typically are distinguished by their average velocities, these velocities may be 
used to roughly classify the asymptotic states as exemplified in Fig. 5. A precise 
classification would require the additional knowledge of the molecule's orientation and 
six bond lengths and hence would be difficult to visualize. 

In addition to the average center of mass velocity the behavior of the chirality 
index Go of the initial configuration is depicted as red line in Fig. 5 along the respective 
cuts. This index correlates with the average velocities at best in a qualitative way. Of 
particular interest here is the behavior of the average velocity in the vicinities of the 
zeroes of the chirality index as a function of the screw angle ip indicated by the blue 
lines. 

As mentioned at the end of Section 2, one has to distinguish between chiral and 
achiral zeroes of the chirality index Gq. Indeed, close to chiral zeroes, none of the 
there existing velocity branches does cross the zero line. But also near the achiral 
zeroes, branches with finite velocities may exist. They may coexist with a branch with 
vanishing average velocity, which, however, may also be missing. The x -B- — x-inversion 
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Figure 4. The center of mass motion of a spiral molecule with an equilibrium 
configuration on cut A specified by h/p = 0.2073, p = 7 and <p — 0.877T asymptotically 
reaches one of the four different states depicted in the panels (a), (b), (c), (d) depending 
on the initial orientation of the molecule. If this orientation is randomly drawn from 
a uniform distribution, then one finds the cases (a) with a probability of 0.19, (b) 
with 0.22, (c) with 0.50, and (d) with 0.09. The ^-component of the center of mass 
is displayed in black. It is shifted by its actual value taken at the time t=3900; the 
y-component, which is either constant, or periodically oscillating is depicted by the 
blue line; from the z-component (green) an average uniform motion is subtracted. 
The configurations mirrored at the y-z-plane starting from the corresponding mirror 
images of the initial orientations of (a)-(d) result in the motion patterns displayed 
in the panels (ma)-(md), which obviously coincide with the mirror images of panels 
(a)-(d). The average velocities (v x ) in ^-direction, see eq. (20) below, in the cases (a) 
and (ma) are ±0.02974; correspondingly one finds (v x ) — ±0.02594 in (b) and (mb), 
(vx) = ±0.00994 in (c) and (mc), and (v x ) = qpO. 03044 in (d) and (md). In all cases 
the mirrored patterns occur with the probabilities of the corresponding original ones 
as specified above. 




Figure 5. Asymptotic average velocities (v x ) were estimated for molecules with 
equilibrium configurations along the four cuts A - D in the <p-h/p parameter plane 
displayed in Fig. 2(a). The initial conditions were randomly oriented according to 
the uniform measure on the three-dimensional rotation group. Panels (a), (b) and (c) 
correspond to the horizontal cuts A, B, C, along h/p — 0.2073 with p=7, h/p = 0.584 
with p = 7, as well as h/p — 2.5162 with p = 3, respectively. Panel (d) represents the 
vertical cut D for <p = 2ir/3, p = 5, and h/p € (0,5). For comparison, the chirality 
index Go defined in eq. (19) is indicated as a red line with the corresponding scale 
at the right-hand side ordinate. The vertical blue lines in panels (a) and (b) mark 
the respective values of the screw angles at which the equilibrium configurations are 
achiral. In all panels the black dots represent the result of the solution of the equations 
of motion (11) with the full Rotne-Prager tensor (6). The average velocities group into 
up to four well defined branches at most parameter values. In panel (a), additionally 
the numerical results obtained from eq. (11) with the Oseen tensor (25) are displayed 
as small green circles. The overall agreement of the two approaches is very good, see 
also Section 4. 



symmetry though is restored on average since there are always pairs of branches with 
opposite velocities which are populated with equal probability. As we shall discuss in 
more detail below, the branches with finite average velocities at achiral zeroes of the 
chirality index represent states of dynamically broken symmetry. 

3.2. Rotational motion 

In the asymptotic state, when the center of mass undergoes a periodically modulated 
uniform motion, the instantaneous rotation axis predominantly points into the re- 
direction, possibly superimposed by small, oscillatory deviations orthogonal to this 




Figure 6. The components of the instantaneous axis of rotation, u>, for the 
configuration with ip = 0.877T, h/p — 0.2073 and p — 7 corresponding to the maximum 
of Go on cut A. The component uj x pointing in the vorticity direction of the shear flow 
is marked in black, the y- and z-components in blue and green, respectively. Panels are 
labeled as in Fig. 4. In all cases the rotation about the vorticity direction is dominant. 
Only in panel (d) a slight tumbling of the axis is noticeable. 



direction. Fig. 6 depicts the components of the instantaneous rotation axis for different 
asymptotic motional states of the same molecule. This orientation of the instantaneous 
rotation axis appears as a rather natural consequence of the considered flow geometry 
with uniform vorticity pointing into the x-direction. A conspicuous tumbling behavior 
of the instantaneous rotation axis was found in a case of dynamically broken symmetry, 
see Fig 12a. The rotation speed undergoes pronounced oscillations in the majority of 
cases. In the asymptotic states in which the instantaneous rotation axis points in the 
x-direction the molecule is typically found to be oriented such that one of its principal 
axes is aligned with the direction of rotation. The transient behavior and the approach 
to the specific asymptotic states of the angles enclosed by the principal axes and the 
instantaneous rotation axis is exemplified in Fig. 7. 

Due to the lack of any symmetry plane in a chiral molecule, there are two possible 
ways to orient a principal axis along the instantaneous axis of rotation. However, 
we found different velocity branches corresponding to opposite orientations only for 
alignments of the rotation axis with the middle principal axis, and only if the skewness 
S u of the molecule along the axis of rotation vanishes within numerical precision. As a 
measure for the skewness of a molecule relative to the vector of angular velocity u> we 
take the third root of the third moment of the atomic positions x n relative to the center 




Figure 7. The time evolution of the cosines of the angles between the instantaneous 
vector of rotation and the eigenvectors of the second moments of positions, cos6*i = 
u> ■ Ui/\u)\, for the configuration with ip = 0.877T, h/p = 0.2073 and p = 7 is depicted 
for different initial orientations. The asymptotic motion of the center of mass as well 
as the motion of the axes of instantaneous rotations are displayed in Figs. 4 and 
6, respectively, with the panels labeled consistently. The cosine of the eigenvector 
corresponding to the largest eigenvalue is red, to the middle one blue and to the 
smallest one green. In panels (a) and (b) the middle principal axis asymptotically 
aligns with the rotation axis, while the initial orientations in panels (c) and (d) lead 
to alignments with the largest and the shortest axis, respectively. 



of mass projected onto the direction of the rotation axis, e u = a?/|u?| , hence, reading 




The maximum number of branches we detected was four. For rj/p = 0.2073, p = 7 and 
screw angles approximately ranging between (p = 0.8tc and tp = 0.9-7T one observes single 
alignments to the small and the large axes and both orientations relative to the middle 
axis. Fig. 8 represents the time evolution of the skewness for a molecule with initial 
orientations belonging to different domains of attraction. 

3.3. Internal motion and dynamical symmetry breaking 

Also in the reference frame that moves with the center of mass and rotates with the 
instantaneous angular velocity, the molecules generally are not at rest but undergo 
internal motions characterized by changing bond lengths and bond angles. These 
distortions then lead to modulations of the magnitudes of the principal axes as 
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Figure 8. For those initial orientations of the configuration with ip = 0.877T 
h/p = 0.2073 and p = 7 which finally approach two different states both rotating 
about the middle axis, see Fig. 7 (a), (b), the third moments of the positions projected 
onto the instantaneous axis of rotation, S M , defined in cq. (21), relax to zero. For other 
initial orientations the skewness S u asymptotically undergoes rotations, see panels (c) 
and (d). 



exemplified in Fig. 9. The resulting periodic relative volume changes of the molecule 
are presented in Fig. 10. 



3.3.1. Dynamical symmetry breaking In molecules with achiral equilibrium configura- 
tions, a deformation caused by shear forces in general will break its mirror symmetry 
and hence render the molecules chiral. This mechanism is at the heart of the effect 
of dynamical symmetry breaking that can be observed at achiral zeroes of the chirality 
index Go- 

For spiral molecules the three-dimensional achiral configurations are found on two 
curves in the <^-/i/p-parameter plane, reading 

ni/2 



(cos(2<£>) — cos(</?)) 



(h/p) 2 = ^ [(cos(3y?) - cos(y?))] 1/2 , 



2tt/3 < if < An/3 , 
tt/2 < if < 3tt/2 . 



(22) 
(23) 



For a derivation see the Appendix. The branch (h/p)i corresponds to the solid line and 
(h/p) 2 to the dashed line in the Fig. 2(b). Fig. 11(a) displays the asymptotic average 
velocities in x-direction and the corresponding average values of the chirality index for 
the achiral configuration at p = 0.72716077T and p = 7 on the curve (h/ p)\. As for all 
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Figure 9. The magnitude of the principal axes of the second position moments for the 
configuration with tp = 0.877T, h/ p = 0.2073 and p — 7. In the four different asymptotic 
states, which are labeled in the same way as in the previous Figures, the magnitudes 
of the principal axes undergo periodic modifications. In the two cases with positive 
asymptotic velocities, (a) and (b), these modifications are almost identical. Noticeable 
differences occur between the states characterized by rotations about the large and the 
short principal axes, displayed in panels (c) and (d), respectively. 



achiral configurations on this curve, it possess five equally long edges. In the present 
case, the sixth edge is shorter than the others. This molecule is initially oriented in 
such a way that the symmetry plane containing both the short edge and the z-axis 
encloses an angle a with the y-axis, see the Fig 11(b). Depending on this angle a, two 
different pairs of asymptotic states with non-vanishing velocities and chiral indices of 
opposite signs are approached with time. The average chirality index was obtained as 
an algebraic mean over a time window of length 1000 in the asymptotic regime. Because 
a rotation about the z-axis by a = n leaves any of these initial configurations invariant 
the resulting pattern of average velocities and chiral indices is periodic with period n. 
The rotational motion and the orientation of the principal axes are presented in Fig. 12 
for two of the chiral attractors of an achiral molecule. 

The average asymptotic velocities on a part of the achiral curve (h/p)i given 
by eq. (22) were determined for p = 5. Similarly as on the other cuts, equilibrium 
configurations were determined from eq. (1) at 2001 equally spaced angular values 
for (p G [0.677r,7r]; subsequently their centers of mass were translated to the origin 
and randomly rotated about the center of mass. With these initial conditions the 
equations of motion were integrated up to t = 8000. The average velocity (v x ) was 
determined according to eq. (20) for to = 7000, N = 1000 and r = 1. Fig. 13 presents 
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Figure 10. The ratio of the actual to the equilibrium volume for different initial 
orientations of the configuration with ip — 0.877T, h/p = 0.2073 and p = 7 in the 
same order as in Figs. 6-9. For the two asymptotic states (a) and (b) moving in the 
positive ^-direction the volume variations are essentially identical but different from 
both states (c) and (d). 




Figure 11. On the panel (a), different chiral-symmetry breaking attractors for the 
achiral equilibrium configuration at tp = 0.72716077T, h/p = 0.584 and p = 7 are 
characterized by their average velocities (v x ) (blue) and average chiral indices (Go) 
(red). The attractors are approached from different orientations of the molecule. The 
angle a specifying the initial orientation is enclosed by the y axis and the molecule's 
symmetry plane containing the short molecular edge as exemplified in the panel (b) 
which presents a view with the z-axis perpendicular to the plane of projection. For 
a = the y-z-plane coincides with this symmetry plane of the molecule. 
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Figure 12. Different modes of the rotational motion of a molecule with cp = 
0.72 7 1 6077T, h/p = 0.584 and p = 7 corresponding to an achiral equilibrium 
configuration on the (/i/p)i-curve. In the upper row the components lj x (black), uj v 
(blue), and tu z (green) of the instantaneous rotation axis are displayed; the lower row 
presents the cosines of angles between the large (red) , middle (blue) and small (green) 
principal axes and the £-axis. The left column corresponds to the asymptotic state 
with larger positive velocity in ir-direction (v x ) — 0.0025 and the right column to the 
asymptotic state with smaller positive velocity in ir-direction, (v x ) — 0.000083. In case 
of the large positive velocity the molecule performs a tumbling but yet periodic motion 
as can be seen both from the oscillatory behavior of the instantaneous rotation axis 

(a) and the alignment of the instantaneous axis with the cc-axis (c). In the asymptotic 
state of slow motion the instantaneous rotation axis exactly points into the x-direction 

(b) . Asymptotically (for times t > 2360) the large principal axis strictly lies in the x- 
direction (d). In a transient period at earlier times the molecular orientation switches 
between perfect ^-alignments with the large and the middle principal axes (d). 



the resulting averaged velocities in x-direction. The emerging pattern is essentially 
symmetric about (v x ) = and consists of different branches and a region with a cloud 
of apparently random points. The deviations from symmetry are insignificant within 
statistical uncertainty. While, on the branches, the center of mass approaches an average 
motion with constant velocity (v x ), superimposed by short oscillations, in the cloud 
region it asymptotically oscillates with a period which is larger than the width of the 
averaging window, see in Fig. 14. Hence, the resulting averages over a window that 
is shorter than a period are randomly distributed within an upper and a lower bound. 
The true averages extended over a full period are zero within numerical precision. 
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Figure 13. Average velocities in the x-direction for 2000 randomly rotated initial 
conditions taken along the achiral curve (h/ p)\ given by eq. (22) are displayed for p = 5 
and ip e [0.67-7T, 0.87r] in steps of Aip = 1.65 • 10~ 4 ir. The equations of motion (11) 
were numerically solved with these initial conditions up to the time t = 8000 and the 
average (v x ) was estimated by means of cq. (20) for t = 7000, t = 1 and TV = 1000. 
Apart from the region ip <G (0.687T, 0.70667r) where the average velocities seemingly are 
random, they form up to four branches. Coming in pairs of opposite average velocities 
they maintain the achiral symmetry at each point along the achiral curve. On the part 
of the achiral curve with p £ [0.8-7T, ir], which is not shown, the lateral average velocity 
forms a single branch with (v x ) = within numerical precision. 
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Figure 14. The center of mass motion on the achiral curve, eq. (22) (a) for ip = 0.77T, 
and (b) p = 0.71092tt, both for p = 5. The long period T w 2730 in (a) leads to 
an under-sampling of the average (v x ) displayed in Fig. 13 and as a consequence to 
the cloud of random points in the screw angle region between 0.667T and 0.70667T. In 
contrast, the case in panel (b) approaches a uniform motion in the positive x direction 
after a short transient. 



As for the case of chiral equilibrium configurations, also for achiral equilibrium 
configurations the number of branches in the present investigation varies from one up 
to four. In the limiting cases of ip = 2/3ti and ip = ir the configurations on the achiral 
curve are planar and only a single branch with (v x ) = exists. Additional branches 
appear for intermediate values of ip always in pairs with positive and negative velocities 
and opposite chiral indices, while the zero-velocity and zero-chirality-index branch may 
become unstable. Branches with opposite velocities and chiral indices always occur 
with equal weight such that the mirror symmetry of the initial achiral configuration is 
restored. For ip = arccos(— 2/3) ~ 0.7327T, corresponding to a regular tetrahedron, four 
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velocity branches cross each other at (v x ) = 0, see Fig. 13. 

The present phenomenon of dynamical symmetry breaking is different from the 
one that was reported in Ref. [17] where molecules with initially achiral geometric 
configurations, but with different bond strengths, were considered. Hence those 
molecules were chiral from the outset with respect to their mechanical properties. 
We emphasize here that the effect of dynamical symmetry breaking occurs for achiral 
molecules with respect to both geometry and bond strengths. 



4. Role of hydrodynamic interactions 

In order to clarify the role of hydrodynamic interactions for the transport of chiral 
objects in a shear flow, we compared the motion of a spiral molecule with screw angle 
ip = 7r/2, relative rise h/p = 0.584 and radius p = 7 in presence of three different vector 
fields describing the action of the fluid on the molecule at different levels of accuracy. In 
the crudest approximation the forces contribute to the velocity in terms of the Stokes 
law for which the mobility tensor simplifies to 

1 

6nr]a 

Leading corrections in the ratios of the atomic radius to intra-molecular distances are 
contained within the Oseen approximation, yielding 



H Stokes U ( 24) 



jyOsccn _ Stokes _|_ ^ 



j _^ r i,j r i,j 



S~//r 



(25) 

<-,j L 

Finally we compared the Stokes and Oseen results with those following from the full 
Rotne-Prager tensor (6). Since we found only minor differences between the outcomes 
of the Oseen and Rotne-Prager treatment, see Fig. 5(a), we only present results of 
the comparison between Stokes and Rotne-Prager. It is interesting to note that the 
rotational motion and the internal motion of the molecule qualitatively are quite similar 
for the Stokes and the Rotne-Prager treatment. Fig. 15 exemplifies the qualitative 
agreement of the two approximations by comparing the alignment of (a) the middle 
principal axis with the instantaneous rotation axis and (b) the re-component of the 
instantaneous rotation axis. 

Also the deformational motion of the molecule is quite similar, see Fig. 16 comparing 
the magnitudes of the principal axes and the relative volume for the two approximations. 
Hence the rotational and internal motion of a molecule is mainly determined by the 
action of the Stokes forces. 

In the present case of uniform shear flow the sole action of the Stokes forces cannot 
lead to a translational motion in the vorticity direction. This rigorously follows from the 
fact that the internal molecular forces exerted by the intra-molecular FENE-Fraenkel 
springs do not contribute to the motion of the center of mass for a uniform shear flow 
within the Stokes approximation, £\ . Hfj okcs Fj = 0, and hence the center of mass 
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Figure 15. A comparison of the rotational motion of a molecule under the sole action 
of the Stokes forces (blue) and the full Rotne-Prager tensor (red) is presented for a 
spiral molecule with ip — 0.87, h/p = 0.2073 and p = 7 with an initial orientation 
that leads to the uppermost velocity branch (see Fig. 5(a)) within the Rotne-Prager 
treatment. The panel (a) displays the approach of the angle between the instantaneous 
rotation axis and the middle principal axis, showing good agreement of the Stokes and 
Rotne-Prager treatment, (b) The oscillations of the x-component of the instantaneous 
rotation axis follow the same pattern with the same amplitude and slightly smaller 
period under the influence of the Rotne-Prager mobility as for the Stokes forces. 



motion is determined by 

X = J2v(x n ). (26) 

n 

Therefore, in a uniform shear flow, as given by eq. (10), the center of mass cannot move 
perpendicularly to the flow direction within the Stokes approximation and the lateral 
drift consequently vanishes. Only the small, non-diagonal terms of the mobility tensor 
implied by both the Oseen and Rotne-Prager treatment lead to a coupling between 
the translational and the rotational motion and, in this way facilitate a center of mass 
motion perpendicular to the flow direction. For more complicated flow fields with non- 
uniform vorticity such a coupling becomes possible and may lead to chiral separation 
even within the Stokes approximation [1]. 

5. Conclusions 

We presented a case study of the dynamics of deformable chiral objects consisting of 
four spheres in a uniform shear flow. These objects, which we have called molecules, 
consist of four spherical atoms. Internal forces due to relative position changes of the 
atoms were modeled by FENE-Fraenkel forces such that the bond lengths may only 
vary within a prescribed range. In the equilibrium configuration the internal forces 
balance each other. Forces exerted by the fluid on the molecule drive the molecule 
out of this equilibrium state. Due to the small size and mass of the atoms, inertia 
terms in the equations of motion could be neglected and a purely overdamped motion 
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Figure 16. For the parameter values detailed in Fig. 15, the hydrodynamic 
interactions contained in the full Rotne-Prager mobility tensor (red) have little 
influence on the magnitudes of the principal axes displayed in panel (a) compared 
to the sole action of Stokes forces (blue). Also the chirality index presented in (b) 
differs in a small shift of the period, which is slightly shorter, and in the amplitude 
which is larger in the case of the Rotne-Prager treatment. 

was studied. The influence of the flowing fluid on the molecules was described by 
means of the Rotne-Prager mobility tensor; this tensor includes the Stokes force and 
hydrodynamic interactions (up to second order in the atom radius to intra-molecular 
distance ratio) between the atoms of a molecule. If only the Stokes forces are considered 
and hydrodynamic interactions are neglected then the molecule moves in a uniform shear 
flow with its center of mass strictly in the direction of the flow field. Moreover it also 
performs a rotational motion about the center of mass with an instantaneous rotation 
axis pointing into the direction of vorticity. However, the hydrodynamic interactions 
cause a coupling between the translational and the rotational motion of the molecule. 
These interactions let the rotating molecule act as a propeller which drives the molecule 
transversely to the flow field in or against the direction of its vorticity. The direction of 
this component of motion is opposite for a molecule and its mirror images, thus leading 
to different transport properties of enantiomers. 

In general, one finds different asymptotic forms of motion for molecules with the 
same equilibrium configuration. These motional states correspond to different attractors 
of the equations of motion. They can be distinguished by the magnitude of their 
average velocities in vorticity direction and the orientation of the molecule relative to 
the instantaneous rotation axis which typically points into the vorticity direction. The 
maximum number of attractors that we identified in the present study was four; this 
corresponds to (i) two alignments of the instantaneous axis of rotation with the large and 
the small principal axes and (ii) two alignments with the middle axis. We documented 
the typical translational, rotational and deformational motion in detail with the example 
of a molecule with a particular equilibrium configuration. The same qualitative motion 
patterns were observed for numerous other equilibrium configurations but not presented 
here. 
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The reason why we studied deformable molecules rather than rigid objects was not 
only to obtain a more realistic modeling but also because the numerical implementation 
of the constraints defining rigid objects in three dimensions is more difficult than 
allowing for deformations. Moreover we found that the flexibility of the molecules 
may lead to dynamical symmetry breaking for molecules with achiral equilibrium 
configurations. We demonstrated that the deformations of achiral molecules may lead to 
chiral configurations with finite chirality index G accompanied by a finite average lateral 
velocity in or opposite to the vorticity direction. If the initial conditions are sampled 
from a symmetric distribution, containing mirror-symmetry related orientations with 
equal weight, then motional states with opposite chirality emerge with equal probability. 
In this way the mirror symmetry of the total state space is recovered. 

The effects of chirality on the transport properties are relatively weak in the present 
case of a uniform shear flow. This finding is also in agreement with the results of Ref. [17] 
for bond-strength induced chirality. The situation though is expected to drastically 
change for non-uniform shear flows. It is plausible to assume that there, locally, still the 
same principal mechanisms are at work: The local vorticity leads to a rotation of the 
molecule which then acts via hydrodynamic interactions as a propeller pushing itself in 
or opposite to the direction of the instantaneous rotation axis depending on its chirality. 
In non-uniform shear flows molecules of different chirality, however, will be moved to 
different regions of the flow-field where they then are advected with different velocities. 
This effect can significantly enhance the influence of the chirality on transport properties 
and hence, may form the basis of an effective, purely physical enantiomer-separation 
mechanisms. 

Before such practical applications can be envisaged, however, the so far neglected 
influences of thermal noise, hydrodynamic interactions between different molecules of 
same and different chirality as well as between molecules and walls confining the fluid 
should all be taken into account and investigated. 
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Appendix A. Achiral configurations 

By the construction of the spiral molecules, the length d n ^ m of the edge connecting the 
nth with the mth atom only depends on the absolute value of the difference between n 
and m, d n ^ m = l(\n — m\). In the general spiral molecule hence has three equal 

edges of length 1(1), two of length 1(2) and one of length /(3) where 




(A.l) 
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One branch of achiral configurations possesses two equilateral and two isosceles 
faces, see Fig 3(c). They have five equally long edges with 1(1) = 1(2). This condition 
determines the relative rise as a function of the screw angle to read 

(h/p)l = ^(cos(2^)-cos(^)) , (A.2) 

where the screw angle is restricted to 2n/3 < p < An/ 3 in order that the right hand 
side stays positive. On this curve the maximal value of the relative rise is reached at 
(p — 7r and there takes the value h/p — 2/\[3. 

Achiral configurations with four identical isosceles faces result from the requirement 
1(1) = 1(3). The relative rise as a function of the screw angle then becomes 

(h/p)l = |(cos(3v?) - cos(v?)) , (A.3) 

with tt/2 < ip < 3n/2. This function vanishes at tp = tt as well as at the borders of the 
indicated interval, and has a maximum at tp ~ 0.6967T, h/p ~ 0.620. 

The two achiral branches cross at the special angle at which cos tp = — 2/3 yielding 
h/p = a/10/27, where all three edge lengths agree with the each other, resulting in a 
regular tetrahedron. 
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